Truncated Lanczos method : application to hole-doped spin ladders 
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The truncated Lanczos method using a variational scheme based on Hilbert space reduction as well 
as a local basis change (Dagotto et at, Phys. Rev. B 58, 12063 (1998)) is re-examined . The energy 
is extrapolated as a power law function of the Hamiltonian variance. This systematic extrapolation 
procedure is tested quantitatively on the two-leg t-J ladder doped with two holes. For this purpose, 
we have carried out calculations of the spin gap and of the pair dispersion up to size 2 x 15. 



The variational approach of the truncated Lanczos 
method (TLM)tira is based first on a proper choice of 
the local basis, and then on Hilbert space reduction and 
exact diagonalizations by the Lanczos algorithm in the 
successive subspaces. This variational technique is a 
good alternative to the density matrix renormalization 
group (DMRG) calculations^ and the pure exact diag- 
onalization technique (ED) for the study of models re- 
lated to high temperature superconductivity. Another 
well known method to study numerically the low energy 
physics of fermionic models on medium size flusters is 
the Quantum Monte Carlo algorithm (QMC)a hindered 
by the famous sign problem. On the other hand, the 
DMRG calculations provide only properties of the lowest 
energy sector and operate on lattice systems with specific 
boundary conditions. The TLM goes beyond these tech- 
nical limitations and allows us to use the ED technique 
on larger fermionic systems. This approach is very use- 
ful e.g. to study momentum-dependence of the observ- 
ables and calculate dynamical properties of Hubbard- like 
models or t-J model on large size systems (compared to 
standart ED techniques) of various geometries. The is- 
sue we shall adress here is the control of the convergence 
of the energy. For this purpose, we will exploit an in- 
teresting behaviour of the energy with the Hamiltonian 
variance, defined as the measure of the fluctuations of 
the energy around its average. We will show that the 
energy extrapolates with a second order power law func- 
tion with the variances, such that the extrapolations to 
the zero-variance limit are accurate estimations of the 
approximate variational energy in a specified sector of 
symmetry. This unbiased extrapolating procedure pro- 
vides more control on the accuracy of the results _ The 
applications to the 2 x L hole-doped spin ladderO (for 
which many exact calculations can be performed up to 
L = 13) show that the error made on the exact val- 
ues are small compared to the degree of approximation 
made. This article begins first with a detailed description 
of the method, followed by a discussion on the efficiency 
of the TLM through applications to the two-hole doped 
t-J model on the 2 x L ladder geometry with periodic 
boundary conditions for the sizes L = 13, which can be 
handled exactly, and L = 15. 



a. Geometry: We will consider closed systems that 
can be formally described like a chain of a finite number 
of interacting finite size sub-systems or clusters. Each 
cluster S 1 is described by a numerically soluble cluster 
Hamiltonian H and two near-neighbor clusters <S' and 
<S J are interacting through the transverse Hamiltonian 
H 1 ^. The cluster Hamiltonian operator Hq acts on the i th - 
cluster Hilbert space £ while the transverse Hamiltonian 
operator H'j_ acts on the direct product £q<E)£ q . Typically, 

contains Heisenberg (magnetic couplings) and kinetic 
interaction terms (one-electron hopping integral). For 
example, in two chain ladders, that we will consider in 
the following, S 1 is simply a rung with two lattice sites. 
The Hamiltonian of the global system is thus a sum of 
two terms H = H + Hj_ with H = i=1 jy- Hq and H^ = 

X)<i j> Hjl , A" being the number of clusters constituting 
the whole system. 

b. Change of basis: We begin with the exact di- 
agonalization of H of all the clusters <S'. The eigen- 
states of the i th -cluster Hamiltonian H are classified 
according to their quantum number (filling, transverse 
and longitudinal momentum, S z spin component). Let 
&h = {l^p >} p=i,...,t , o W1 th 'Po = dim£ denote the lo- 
cal basis set of the eigenstates of H in £q. The direct 
product £ = £q ®£q <g> ... ®£(f of the A" cluster Hilbert 
spaces £g defines the Hilbert space of the whole system. 
The initial basis set Bo of £q in a given sector of symmetry 
is thus defined by the following set of A/"-tensorial product 

vectors ^\& q >= ®i=i,jv l^pi > } reduced to the states 
consistent with the quantum numbers implemented for 
the global system (in practice Qo — card So -C V$). 
Typically, use of translational symmetry in the longitu- 
dinal direction provides a significant reduction of Bq. We 
notice that Ho is naturally diagonal in Bo- Its eigenvalues 

are just the sum ho,{Pi,-,PAf} = h o,Pi + n o, P2 + - + h ^ 
where the set {pi, ...,paa} is labelling a particular set of 
cluster Hamiltonian's eigenstates (h is the eigenvalue 
of H associated to \<j? >). 

c. Expansion and truncation of the Hilbert spaces: 
We have represented on Fig. |l| a picture of the itera- 
tive scheme. Although Bq can be formally constructed it 
would be impossible to diagonalize exactly the Hamilto- 
nian in this basis for sizes of interest, so that a truncation 
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procedure is necessary. Keeping a small identical subset where I„' is the set containing the indexes of the ba- 
sis states |$ q > that belong to the reduced basis B' n . 



of the states in all the local basis B Q would be in fact a 
quite poor approximation and we shall use a better pro- 
cedure. The first step consists in choosing a subset B' 
of Bo . A simple choice is to restrict ourselves to the few 
low energy states of B . The exact diagonalization by 
the Lanczos algorithm of the Hamiltonian H = Ho + H_l 
in B' gives the energy e' in a fixed sector of symme- 
try at the zero th -order of the iterative expansion. The 
zero th -order approximation of the ground state is then a 
normalized linear combination of the vectors in B' . At 
the first order, the expanded variational basis B\ is cal- 
culated by applying the transverse part Hi on the states 
in B' (card Si » card£> ). Like previously, we calculate 
the first order of the variational energy ei by an exact 
diagonalization of H in the Hilbert space £\ generated by 
B\ (of course ei < e' ). Then the best states (see later) 
in B\ which constitute the new reduced basis B[ are se- 
lected. The reduction of the Hilbert space £ i ends the 
first iteration. The rest of the calculation has the same 
loop-structure. Let us briefly recall the scheme of the 
(n — l) th iteration (Fig. [[]). It consists in two steps : (i) 
Expanding B' n _ 1 by applying H_l. Let B n denote the new 
expanded basis (cardZ?„ 3> cardZ?'^), (ii) Truncating 
the generated basis B n , which give the reduced basis B' n . 
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The DMRG method also use such approximations. In 
order to select the best states we define a parameter 
e £ [0, 1], such that a state |$ q > belongs to the reduced 
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basis if the associated weight obeys to 
The actual value of e at each iterative step can be de- 
termined by three different methods : (i) by specifying 
the minimum value of the overlap O = < ^ n \T^n > 
~ SqGi ' II a ^ II 2 ' fixing the fraction r\ of the 

states of B n kept in the reduced basis set B' n and (Hi) 
by fixing an expansion parameter A > 1 such that 
dim£' n — A ■ dimS^—v In the following, procedure (Hi) 
is used (for n > 2) since it gives more control on the 
size growth of the Hamiltonian matrix and on the CPU 
time of the calculation. The first reduced space £[ is 
obtained by using the criterium of minimum overlap, ie 
< *i|T*i >= 0\ > O min . Note that for moderate ex- 
pansion parameter A (see table |), a rather large fraction 
of the weight is preserved at each step. 



n th loop 



A = 2.27 
dim £' n dim £ n 



2 048 33 804 460 
0.9990 10 278 21 504 
0.9967 23 299 83 679 
0.9959 52 939 169 045 

0.9961 119 996 327 084 



O 



X = 1. 
dim£ ' 
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dim£ n 



0.9991 
0.9927 
0.9869 
0.9812 
0.9761 
0.9772 
0.9826 



2048 
10 278 
15 476 
23 320 
35 099 
52 846 
79 599 
119 967 



33 804 460 
21 504 
83 679 
126 031 
174 257 
245 066 
347 059 
472 799 



TABLE I. Evolution of the Hilbert space dimensions with 
the truncation criteria A = 2.27 and A = 1.51, for the 2 x 13 
two-hole doped spin ladder in the (S z = 0, K = (0, 0)) sector. 
We reported also in the first column the overlap O (defined in 
the text). The basis B' is built with the 2048 lowest energy 
eigenstates of Ho in Bo (with two holes on the same rung 
and particles on the other sites) compatible with the global 
quantum numbers. 



FIG. 1. Schematic picture of the successive Hilbert space 
truncation and expansion procedures. The operator T acts 
on the Hilbert spaces £ n to obtain the reduced spaces £' n . We 
described precisely in the text the action of the non-unitary 
operator T. 

The normalized variational eigenvector > at order 
n of the truncation procedure is a linear combination of 
the Q n = dim £ n vectors of the subset B n of Bo ■ In order 
to reduce the dimension of the Hilbert space £ n ( n > 1 ) 
we use a systematic procedure to select the best states in 
l^n >. Let T denotes the non- unitary truncation opera- 
tor defined by the following Hilbert spaces morphism : 



T : £ n 

|*n > 



£' 

|T*„> = 



E 

q ei„' 



(1) 



d. Extrapolation procedure: As it will be clear in the 
following, it is very useful to define the Hamiltonian vari- 
ance a n at the n th iteration : 
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{<<|h 2 |*: 1 >-<*:jh|*; 1 > 2 } (2) 



where N is the number of sites on the lattice and the 
new wave function \^' n > is the normalized eigenstate of 
H in the truncated basis B' n with the lowest eigenenergy. 
The calculation of a n does not need any approximation 
and can be exactly evaluated. < \ H 2 | *S?' n > is cal- 
culated by expressing the vector H > in B n +i and 
evaluating its norm, while < ^' n \ H \^f' n > 2 is easily ob- 
tained by the exact diagonalization of H in B' n . o~ n mea- 
sures the quadratic deviations of the energy from its aver- 
aged value at the n th iteration. We have tried four kinds 
of extrapolation of the energy versus the Hamiltonian 
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variance : (i) e n vscr n , (ii) e' n vsa ni (Hi) e n +ivsa n and 
(iv) e' n+1 vsa n where e„ =< *„| H |*„ >, with |* n > 
the normalized eigenwave function associated to the low- 
est eigenenergy of H in £>„, and e' n — < ^f' n \ H \ W n >. 
The zero-variance limits of the extrapolation curves are 
approximations of the exact energy. Note that each value 
of the variance corresponds in fact to a distribution of en- 
ergies since two variational states of same energies might 
have different variances. Although we did not explicitely 
evaluate this distribution of energy, this uncertainty is 
directly connected to the (estimated) difference between 
c' n and e„ for the same variance. Hence, this gives a 
procedure to take into account this uncertainty in en- 
ergy while extrapolating our results and we shall discuss 
quantitatively the precision of this technique through ap- 
plications to the t-J model defined on the two-leg ladder 
geometry. 
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FIG. 2. Energy of a 2 x 11 (a), 2 x 13 (b) and 2 x 15 (c) 
ladder with two holes extrapolated as a quadratic function of 
the Hamiltonian variance compared to ED results when avail- 
able. We made the calculations for 6 iterations with a final 
Hilbert size fixed to 120 000 symmetric states in the sector 
[S z = 0, K = (0, 0)) for the (b) and (c) cases, while in (a) we 
plotted the results for 8 iterations with a final Hilbert space 
quite equal to the complete Hilbert space of the 2 x 11 lad- 
der with two holes (1 940 064 states without the spin inversion 
symmetry) in the sectors (S z = 0, K = (K x ,0)) with K x = 
and K x ~ 7T. 

e. Applications to the 2 x L hole-doped ladder: In 
this section, we apply the TLM to the t-J model on the 
2x L Heisenberg spin ladder with two holes and periodic 
boundary conditions. The Hamiltonian is defined by 

H = J y^(Sj,A • S i+ i,A - -n^xfii+i,),) (3) 
+ J y^(S^i ■ S h2 - jn hl n iy2 ) 



i.s 



(4,!, 



+ h.c.) 



Here c\ < creates an electron of spin s on site i of 
leg A = 1 or 2, S i)X = (4a,s cr ss /C i A s ,)/ 2 and «t,A = 
S s c| A s c i x s- We have taken both the nearest-neighbor 
leg and rung one-electron hopping matrix elements equal 
to t — — 1. The exchange interaction J is taken as 
isotropic between nearest-neighbor leg and rung sites as 
well and throughout J/t = 0.5. Extrapolations of the 
energy versus the Hamiltonian variance are shown on 
Fig. ^ for the 2 x 13 and 2 x 15 two-hole doped ladder 
in the sector (S z — 0, K = (0,0)). These extrapolations 
are mainly linear, with a quadratic correction at low vari- 
ance. The energies corresponding to the zero-variance 
limit on these curves should be compared to the exact 
ground-state (GS) ones Eqs (for sizes 2 x 11 and 2 x 13 
see RefE). For the 2 x 13 ladder, we noticed first that 
there is no significative difference between the extrapola- 
tions corresponding either to 5,6,7 and 8 iterations, but 
also for those associated to final Hilbert space sizes of 
120 000 and 200 000 states. The stability of the extrapo- 
lations with variations of these parameters is an indica- 
tion of a good precision on the approximate energy. The 
two best approximations of the variational energy come 
from the extrapolations of e„+i and e' n+1 vs er„. We min- 
imize the error by averaging the two values of the energy 
obtained by these extrapolations. For the 2 x 13 sys- 
tem, this procedure leads to an error on the energy less 
than 1.85% of the energy gap between the first excited 
level and the fundamental (0.005 compared to ~ 0.27), 
while we took into account only 0.37% of the states in 
the Hilbert space of the whole system. This means that 
only a negligible fraction of the basis states in the funda- 
mental wave function have a non-negligible contribution 
to the total weight. The calculations on the 2 x 15 lad- 
der with the same parameter A generate an averaged ex- 
trapolation with more uncertainty (~ 0.2) than the one 
observed for smaller system sizes. This is mainly due to 
the decreasing of the overlap of the truncated states at 
each step of the calculation, so that we should increase 
the dimension of the final Hilbert size in order to com- 
pensate this loss in precision. Meanwhile, the evolution 
of the error bars give us an information on the accuracy 
of the extrapolated values of the variational energy, in 
such a way that we are able to adjust the initial para- 
maters to reduce the error margin. We have attempted a 
systematic extrapolation procedure within this method. 
As a test, we will complete calculations of the spin gap 
and of the pair dispersion of the low lying S z = pair 
state for two-hole doped 2 x L ladders up to size L = 15. 
As in a preceeding work onpthe evolution of the spin gap 
upon doping a 2-leg ladderQ we defined the spin gap as 
A s - E Q (n h = 2, S = 1) - E {n h = 2, S = 0). In this 
paper, we used the ED Lanczos method on 2 x L ladders 
of maximum size L — 13. With a less important nu- 
merical effort, we expect improving the precision of the 
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TLM by comparing the spin gap obtained with the TLM 
for L — 13 with the exact calculations. Our data follow 
in general a power law behaviour (Fig. ||). Note that in 
the sector (S z — 0, K = (K x ,0)) the extrapolations are 
dominated by the linear behaviour for low momentum 
while they are rather quadratic for higher momentum 
in the same sector and for (S z = 1, K = (tt,tt)). Thus, 
the spin gap appears logically to be extrapolated by the 
polynomial function As = ao + a% <r n + a 2 The aver- 
aged values in the zero-variance limit are As = 0.14 for 
L = 13 and As = 0.10 for size L = 15, while the exact 
value is As° xact = 0.19 for L = 13. The lost in accuracy 
of As for L = 15 is connected to the increased uncer- 
tainty of the average of the extrapolated values due to 
the decrease of the proportion of states kept in the last 
truncated Hilbert space. 
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FIG. 3. Extrapolations of Ae = (e„ + e' n )/2 — where 
eoo is the zero-variance limit of the averaged GS energy for 
the 2 x 13 and 2 x 15 ladders with two holes in the sectors 
(S z = 0, K = (0, 0)) and (& = 1, K = (n, tt)). We made the 
extrapolation of the energy calculated for 6 iterations with 
120 000 symmetric states in the last truncated Hilbert space 
(0.37% and 0.0021% of the whole Hilbert space of the 2 x 13 
and 2 x 15 ladders). As is obtained as the zero- variance 
extrapolation of the upper curves and compared to ED result 
for 2 x 13. 

Lastly we apply this systematic extrapolation pro- 
cedure to study the pair dispersion in the sector 
(S z = 0, K y = 0). The results obtained within the TLM 
for L = 11, L = 13 and L = 15 are shown on Fig. |J. For 
the 2x11 and 2 x 13 ladders with two holes, we plot the 
variational energies e„ obtained at the n = 7 th iteration 
for final truncated Hilbert space sizes fixed respectively 
to 120 000 and 1 300 000 states, which correspond approx- 
imatively to the same proportion of the complete Hilbert 
space sizes (6.5% and 4% respectively). In the case of 
the 2 x 15 ladder, we kept 120 000 states (~ 0.0021% of 



the whole space) and we made 6 iterations. Let us recall 
that the error bars correspond to the difference between 
the extrapolations to the zero- variance limit of e ra +i and 
c' l+1 in function of a n . The loss of accuracy on the aver- 
aged value of the extrapolated energies is obvious when 
the ladder size is increased for the same number of states 
in the last truncated Hilbert space, and we notice also 
that the error on the extrapolated energies for the 2 x 15 
ladder is also increasing from 0.2 to 0.6 for K x > ir/2. 
For small sizes, the extrapolated pair dispersion coin- 
cides with the exact one in the low momentum sector 
(K x < 7r/2). For high momentum (K x > ir/2), we no- 
tice on the three figures a separation between the pair 
dispersion obtained with the variational energies which 
remains close to the exact one and the average of the 
extrapolated energies. This fact is probably related to 
the proximity of the continuum of the excitated states in 
high momentum sectors. The dispersion of energies for 
each value of the variance is increasing, so that our punc- 
tual extrapolations become less accurate than in the low 
momentum sector. Clearly this point needs to be more 
precisely elucidated. 
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FIG. 4. Pair dispersion calculated with the truncated Lanc- 
zos method on the 2 x 11, 2 x 13 and 2 x 15 ladders with two 
holes in the sector (S z = 0, K y — 0) (see the text for details). 

To conclude, we have described and tested a trunca- 
tion procedure which can be applied to a large class of 
strongly correlated systems. This is a tool to construct 
variational approximations of the exact ground state by 
only keeping a small fraction of the whole Hilbert space. 
Here, we have focussed on the reliability of the method 
to produce accurate energies, an issue not addressed in 
details in the litterature so far. We have discovered an in- 
teresting behavior of the variational energy as a function 
of the variance of its corresponding state which enables to 
quantitatively estimate error bars. The method has been 
tested in the case of a two-hole doped two-leg ladders and 
we have found that the accuracy depends strongly on de- 
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tails of the excitation spectrum (presence of a gap, of a 
continuum...). In particular, the accuracy of the two-hole 
pair dispersion is better for momenta where the ground- 
state is well separated from the continuum. Although, 
the method is certainly not as accurate as e.g. the DMRG 
technique, we believe it can offer a complementary anal- 
ysis in case where the DMRG is less efficient (e.g. for 
systems with periodic boundaries) while enabling to use 
symmetry analysis as in standard exact diagonalisation 
methods. Systems where part of the charge or spin fluc- 
tuations are selectively suppressed by external potential 
(e.g. stripe potential, etc..) are best candidates for this 
method. 

We thank IDRIS (Orsay) for allocation of CPU time 
on the NEC SX5. 
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